Association of heteronuclear molecules in a harmonic oscillator well 
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We describe the collisional interaction between two different atoms that are trapped in a harmonic 
potential. The atoms are exposed to a magnetic field, which is modulated in the vicinity of an s-wave 
Feshbach resonance, and we study the formation of molecular bound states and excited states of 
the trapped system with non-trivial angular correlations. 
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I. INTRODUCTION 

Among the variety of processes and phenomena that 
have been studied with degenerate quantum gasses, the 
conversion of an atomic Bose P, |3, S iJ] or Fermi- 
gas d, H, Q of atoms into a gas consisting of diatomic 
molecules is one of the most fascinating, because the en- 
tire system undergoes a dramatic change of state and be- 
cause it provides a practical way to produce a quantum 
degenerate state of a molecular species, which may not be 
reached easily in any other way. This process, which can 
be driven both by photo-association, by sweeping a B- 
field across a Feshbach collisional resonance in the system 
and by RF association from another hyperfine state, has 
been studied both in larger trapped samples where the 
collective many-body state changes character, and in the 
Mott insulating phase in optical lattices 0, [1, B [13, [HI , 
where the process involves only two particles and can be 
understood microscopically. 

With the association of different atomic species one has 
the possibility to form polar molecules, which offer in- 
terestin g in teraction dynamics both in large condensates 
[13, [3 13 [3 J and in Mott insulators, where the long 
range dipole-dipole interaction between molecules, e.g., 
makes them an interesting candidate for quantum simula- 
tors and quantum computers [TgI . [l7j . In the present the- 
oretical paper we study the dynamics of a pair of different 
atoms in the vicinity of an s-wave collisional resonance. 
We treat the case of atomic systems prepared in an op- 
tical lattice potential with precisely one atom pair per 
site, and we treat the case of a deep lattice well approxi- 
mated by a harmonic potential. This system has recently 
been implemented and studied [7| , and the production of 
molecules by RF association from another hyperfine state 
was demonstrated. We shall identify the eigenstates of 
the problem for different values of the B-field and solve 
the dynamical equations when the field is modulated. In 
particular we will show that oscillations of the field with a 
frequency which is resonant with the Bohr frequency be- 
tween discrete atomic-like and molecular-like levels in the 
trapping potential offers unique opportunities for control- 
lable transfer of the atoms to specific excited and molec- 
ular bound states in the trap. 



In Sec II, we outline the basic theory, which is compli- 
cated here by the fact that the problem with different 
atoms does not separate in center-of-mass and relative 
coordinates. In Sec. Ill, we treat the atomic interaction 
by the conventional pseudo-potential, and we show that 
all necessary matrix elements of the Hamiltonian can be 
determined analytically, and the eigenstates and energies 
can be found by diagonalization of not too large matri- 
ces. The eigenstates have definite total angular momen- 
tum, but the non-separability implies that the states are 
not eigenstates of the separate center-of-mass and rela- 
tive angular momentum. In Sec. IV we generalize our 
earlier dynamical calculations [l^l to the heteronuclear 
case. Using a harmonic oscillation of the magnetic field 
we demonstrate association of heteronuclear molecules 
with 100 % efficiency and show that even though we start 
from spherically symmetric states and have only s-wave 
interaction and symmetric external potentials, it is possi- 
ble to drive the system into states with non-zero relative 
angular momentum. Sec. V concludes the paper. 



II. NON-SEPARABLE MOTION OF 
DIFFERENT ATOMS IN A HARMONIC 
POTENTIAL 

If the potential depth in an isotropic 3-dimensional 
optical lattice, V{r) — Vq X]i=x y 2 ^i'^^ (^^0, is signif- 
icantly larger than the tunnelling energy between lat- 
tice sites, atoms placed in the lattice may be prepared 
in a Mott-insulating state with a well defined number 
of atoms in each well. Provided that the lattice is also 
deeper than the recoil energy h^k"^ /2ra, the atoms expe- 
rience around each potential minimum a harmonic oscil- 
lator potential 



V{v) ~ Vot 



i—x^y,z 



Vok^r^ 
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(1) 



where w = •y/2Vbfc'^/m, k is the optical wave number and 
m is the mass of each atom. 

We shall treat the dynamics of two atoms, with masses 
nil and m2 and position coordinates ri and r2, from 
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which we form the center-of-mass (CM) and relative co- 
ordinates 

TOiri + TO2r2 

R= ■ r = ri-r2. (2) 

mi + 1712 

The Hamiltonian with an atomic interaction potential 
Vint(r) reads 

h"^ 1 1 

H = - — - —Vl + -mu;lrl + -mulrl + M„t(r), 
2mi 2m2 z 2 

(3) 

and it does not separate in center-of-mass and and rela- 
tive coordinates, but it can be rewritten 



with 



if = i/cm(R) + £rrci(r)-HC7R-r, 



2M 



(4) 



(5) 



1 



i/rcl - -^^l + -M^V^ + Mnt(r). (6) 

In ([5]), we have introduced the total and reduced 
masses M = mi + m2, = mim2/M. 

The frequencies and w and the coupling coefhcient 
C have non-trivial values determined by equating the po- 
tential terms in ^ and 



1 + /3 



1 + a ' LUi 



+ P/a C 



1 



miUJt 



1 



(7) 



where we define the mass ratio a — 1122/ mi and the ratio 
of the harmonic oscillator depths (3 — {m2U!2) / i^i^i) ■ 

Note that in the case W2 = (that is, a = /?) we 
have ^ = uj = UJl = lu2 and C = so the center-of- 
mass and relative motion separate exactly, independently 
of the masses. In the general case we obtain a simple 
CR • r coupling of the two degrees of freedom, which 
suggests to diagonalize this coupling in the factored basis 
of eigenstates of the separate problems. 

Let us take ®^Rb and ^Li as a specific example 
with a rather large mass ratio. If we take (3 = 
(mLiWLi)/('^Rb'j-'p{j) — 0.4 (which means that luu = 
2.41 X WRb) we have ft = 1.14 x Ci^Rb, uJ = 2.34 x oj^h, and 
C = —0.877 from ([7]). Taking the harmonic oscilla- 
tor length scale as an estimate of the magnitude of the 
dipole matrix elements, we can get the following order 
of magnitude estimate of the energy shift of the lowest 
energy states caused by the coupling: 




\C\ 



0.3 



(8) 



so the coupling is expected to be significant but the prod- 
uct basis is still a good starting point for our analysis. 



We shall assume that the interaction potential is cen- 
tral, i.e., it is independent of the relative angular orien- 
tation of the atoms. This implies that the iJcm and iJrci 
commute with the angular momentum operators, 



L = R X P 



Z = r X p 



(9) 



where P and p are the conjugated momentum operators 
of R and r respectively and we can choose as basis for 
our description the product states 



NLMnlm 



(R,r) = $ArLA/(R)V'nim(r) 



(10) 



where N and n are principal quantum numbers and L, M 
and /, m are the angular momentum quantum numbers of 
the CM and relative motion respectively. It is interesting 
to note that the total angular momentum of the particles 
J = li + I2 can also be written as 



J = L + L 



(11) 



Since both Hcm, H^d and R • r are rotationally invari- 
ant, the total angular momentum commutes with the 
full Hamiltonian, and our Hilbert space separates into 
independent subspaces belonging to different values of 
J and mj which makes it natural to switch from the 
\NLMnlm) basis to the coupled \Jm,jN Lnl) basis: 

^iTnl = Y.^LMlm\Jmj)<^NLM{'R)^nUn{T^) (12) 

Mm 

where {LMlm\Jmj) denote the Clebsch Gordan coeffi- 
cients. 

To diagonalize the R • r coupling term in the basis (fT2l) 
we need the matrix elements: 



J 771 J 



{L'M'l'm' \Jmj) {LMlm\Jmj) 



M' in' Kl7n 



(r)|R • r|$ArLAf(R)<(5„im(r)) 

(13) 

It should be emphasized that R-r does not commute with 
and F so L and I are not good quantum numbers when 
the coupling term is included into the Hamiltonian. This 
is in contrast to, e.g., the L-S (spin-orbit) coupling in the 
hydrogen atom which commutes with Lp' and leaving 
L and 5* as good quantum numbers. 



III. EIGENSTATES OF RELATIVE MOTION 

A. Pseudopotential description 

At sufficiently low temperatures atoms interact mainly 
through s-wave scattering since the centrifugal barrier for 
higher partial waves becomes much larger than the ther- 
mal energies. The interaction depends to a large extent 
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only on a single parameter, the s-wave scattering length 
a^c, and the real physical interaction can be modelled by 
the following pseudopotential [l^ [1^ 



^ or 



(14) 



where the action of the operator {d/dr)r on a wave func- 
tion tp is to be understood as {d/dr){rip). The pseudopo- 
tential reproduces the correct wave function in the entire 
range outside the range of the physical interaction. 

Writing u{r) — ri(;{r), the pseudopotential implies 
the boundary condition plj |. 



u'jO) _ -1 

m(0) flsc 



instead of the usual u(0) = in a regular potential. 

For interatomic potentials containing a long range van 
der Waals {V{r) = —Ce/r^) attraction, the validity of the 
pseudopotential description for trapped atoms in a har- 
monic oscillator is determined by the ratio Pe /oho where 
(3q — {2^Ce/h'^y/'^ is the characteristic potential length 
scale and aho = \/WJjio is the length scale of the har- 
monic confining potential [H, [2^ . Values of Ce for both 
homonuclear and heteronuclear alkali metal atom pairs 
can be found in [23|. In the case of an optical lattice 
well with a trapping frequency of = 27r x 200 kHz we 
have Cho = 1-2 x 10"^ oq and since /Jg — 85 ao for *^Rb-^Li 
[2^ we expect the pseudopotential approximation to be 
reasonably good even for a very tightly confining lattice. 

The real interaction potential accommodates several 
bound states, but here we are only concerned with the 
very loosely bound state which can be accessed in Fes- 
hbach resonance experiments and this bound state is 
accounted for by the pseudopotential. The parameter 
which is used to tune the scattering length in experi- 
ments is the magnitude of the magnetic field B. Near a 
Feshbach resonance the scattering length is given by [2^ 



1 - 



B- Br 



= flbg 1 - 



B-Br 



(16) 



so if everything is expressed in terms of the dimensionless 
magnetic field [B — Bo)/A the only tunable parameters 
in the pseudopotential model are the harmonic oscillator 
frequency uj and the background scattering length abg 
(notice, however, that A might be negative). 

Calculations of spectra in models with more accurate 
potentials can be found in [g^, H^, [13] , and pseudopoten- 
tial models with an energy dependent scattering length 
have been used to improve the accuracy [12, [11| . While 
the exact spectra and matrix elements might be slightly 
different due to the limitations of our approximation, the 
pseudopotential captures the essential physics quite well 
and it gives good quantitative agreement in most cases. 



B. Wave functions of relative motion 

In this section we obtain solutions to the Schrodinger 
equation of the relative motion with the Hamiltonian 
© and the pseudopotential These solutions, i.e., 

the energy spectrum and the wave functions have been 
provided previously [U [1^ , but we give them here for 
completeness, as we shall use their precise form to sub- 
sequently evaluate the matrix elements of the coupling 
term CR • r. 

Let us first, for reference provide the spectrum 



E = nLu[2n + l 



(15) and the associated eigenfunctions 



V'nim(l") 



2n! 



r(n + / + 3/2) 



(17) 



(18) 



(19) 



for the three-dimensional isotropic harmonic oscillator. 
L" are generalized Laguerre polynomials, r is given in 
units of the harmonic oscillator length yjh/{^u!), and 
Yim are the spherical harmonic wave functions. 

The pseudopotential only alters the solutions with 
I — [28]. In the region r > the radial Schrodinger 
equation, expressed in units of the harmonic oscillator 
length scale Oho = \J^I [i-^ and energy scale E = hue, 
reads 

15^" 1 2 
-2^ + r" = - (20) 

or 

According to [29| 19.3.7 and 19.3.8, the solutions 
to this equation are the parabolic cylinder functions 
Dc-i/2{V2r) and V(—e, V2r) with the following asymp- 
totic behaviour for large r ([2^ 19.8.1 and 19.8.2): 



A-i/2(V2r 
V{-e,V2r) z 



2^.V2 



/2(V2r)'-i/2 



(22) 
(23) 



The diverging V function must be discarded. The D 
function satisfies ([li] 19.3.5 and 19.3.7) 



i5e-l/2(0) 



2-^/2+i/4r(-e/2 + 3/4) 



2-^/2-i/4r(-e/2 + l/4) 



(24) 



(25) 



Normally we would only accept the regular solutions with 
m(0) = 0, which implies that the F function in ([M]) is 
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infinite, and hence that its argument is a non-positive 
integer. Thus, e = 2n + 3/2 where n — 0,1,2,... in 
agreement with (fT7| . The pseudopotential 111]), how- 
ever, imposes the boundary condition (jl5p . and the dis- 
crete energy spectrum is given by the solutions to the 
foUowing transcendental equation: 

V^^Ui/2(0) ^ -2r(-£/2 + 3/4) ^ -Oho 
De-i/2(0) r(-e/2+l/4) a,c 

It is convenient to associate to each eigenenergy an effec- 
tive harmonic oscillator quantum number v defined by 



(26) 



such that u{r) cx D2v+i{^/2r) and 

Mr) - ^Yoo = ^2-^-^/'D,.+,{V2r)Yoo (27) 
r r 

where 



' 2T{-iy-l/2)r{~i^) 



(28) 



is a normalization constant and ip is the digamma func- 
tion (see Appendix). The wave function generally di- 
verges as 1/r for small r since u(0) is generally finite 
according to The expression ([?f|) can also be rewrit- 
ten in terms of the confluent hypergeometric U function 
(see I2p 13.6.36 and 13.1.29) to obtain the forms cast in 



C. Dipole matrix elements and non-separable 
eigenstates 

In a given Jmj subspace the matrix representation 
of the Hamiltonian ^ in the \NLvl) basis, abbreviated 
by state indices p, q, is the sum of diagonal contributions 
from i/cM and H-cci and an off-diagonal contribution from 
the coupling term CR • r (see Eq. [T3| : 

Hpq = Ml{2Np + Lp + ?,/2)5pq + huj{2vp + Ip + i/2)5pq 
+ C ^ {LpMplpmp\Jmj){LqMglgmq\Jmj) 

MpmpMqniq 



X ($A,^L^M,(R)|R|$A',L,Af,(R))R 
■ {^^plprap{Y)\v\ipu,l,-m,{Y))v 



(29) 



where ( )r and ( )r denote integration with respect to 
R, respectively r and 5pq is the Kronecker delta. 

The dipole integrals separate into radial and angular 
components. The radial dipole integral for harmonic os- 
cillator wave functions is given by the following simple 
formula, see Appendix: 



RN'(L+i){r) r RNLir)r'^dr 



= ^/N'TTT3/2 6n',n - Vn 6n',n-i (30) 



The parabolic cylinder wave functions have Z = so they 
couple only to the / = 1 harmonic oscillator states. In the 
Appendix we denote the radial parts of the wave func- 
tions by Ri, and and we derive 



R,ir)rRni{r)r^dT 



2nl 



r(n + 5/2) 16 



k=0 



-k ) k\ T{k + 2-vy 



(31) 



We now have analytical expressions for the full Hamil- 
tonian matrix. Consider the case where the system is 
prepared, and therefore remains, in the block of J = 
0, mj — states. Such states involve CM and relative 
states with L = I and M = —m and using the Clebsch- 
Gordan coefflcient 

{LMlm\J = TOj = 0) = 6M{-m) i-l)'^~^'/V2L 
we have 



V'^"^„i(R,r) 

= E (-l)''-''<fiVLAf(R)v'„L(-M)(r). 



1 



V2L 



M=-L 



(32) 



When we diagonalize the Hamiltonian in the basis ([32]) 
we get the J = spectrum shown in Fig. [TJ Note that wc 
plot the energies as functions of scattering length, which 
in turn must be described as a function of the magnetic 
field in an experiment. 

As the eigenfunctions ^' of the complete Hamiltonian 
are expanded on the basis (|32p of angular momentum 
eigenstates 



(33) 



NLr, 



the expansion coefficients provide, through |cArL„p, the 
CM and relative angular momentum content of the eigen- 
states. We observe that with a central interaction po- 
tential and spherically symmetric single-particle confin- 
ing potentials, it is possible to generate states with non- 
vanishing relative, and center-of-mass, angular momenta. 
In Tab. |T]and in Tab. |n]the population of different angu- 
lar momentum components are given for the lowest eigen- 
states of the system for Osc = and for Csc/aho = 0.36. 
For instance the state with an energy of 3.64 huj at age = 
contains 63.2 % of |010) with L = I = 1, 13.6 % of |100) 
with L = I = 0, and 23.1 % of other product states. 
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Figure 1: (Color online) J = spectrum of a Rb and a Li atom 
in a harmonic well in the pseudopotential approximation. The 
energy scale huj of the relative motion harmonic oscillator is 
chosen as the unit of energy. /3 = (mLiti'Li)/("iH,b'^Rb) ~ 0.4. 
The spectrum of the full Hamiltonian including the coupling 
between CM and relative dynamics is shown as the connected, 
blue lines. The spectrum we get before including the coupling 
term is shown as red, dashed lines. We have included states 
up to A'^ = 6, n = 6 and L = I = 12 which gives a total of 637 
product basis states. Note the equidistant CM vibrational 
states of molecular nature when a^c > 0. 



Energy 


Largest content 
\NLn) |civL„P (%) 


Second largest content 
liVLn) Icivinp (%) 


2.18 


|000) 96.2 


|010) 3.5 


3.04 


|100) 79.1 


|010) 11.2 


3.64 


|010) 63.2 


|100) 13.6 


3.89 


|200) 64.0 


110) 16.7 


4.24 


|001) 80.8 


|010) 12.6 


4.49 


|110) 33.5 


|200) 22.8 


4.74 


|300) 51.1 


|210) 18.3 



Table L Angular momentum basis content at a^c = of the 
eigenstates of the full Hamiltonian for ^'^Rb and "Li in a har- 
monic potential with (3 — (mLitt'ii)/(?nRbCij^(,) = 0.4. 



IV. RESONANT DYNAMICS NEAR A 
FESHBACH RESONANCE 

A recent experiment Q has demonstrated the forma- 
tion of heteronuclear molecules by RF association from 
another hyperfine state near a Feshbach resonance in the 
87Rb-40K system. 

Extending our previous work [isj for identical sys- 
tems we will now address the dynamics beyond the adi- 
abatic approximation, where a modulation of the field 
causes transitions between the adiabatic eigenstates. As 
a specific apphcation we will discuss resonant transitions 
driven by an oscillating magnetic fields The idea, pro- 
posed and used by Thompson et al. I30| and Greiner et 
al. [m and theoretically analyzed in [islls^. is to apply 



Energy 


Largest content 

\N Ln) \CNLn\ (%j 


Second largest content 

\N Ln) \CNLn\ (70j 
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110; 0.002 


11/1 

-1. 14 


IOAA\ AO AAA 
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1 01 A\ A AAQ 

|ziu; U.UUO 


A 1 

-O.lD 


IOAA\ AA AOK 

|300) yy.yso 
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310; 0.004 


A O 1 

O.ol 


1 A AA\ AA AO 1 

|4U0) y9.981 


/I 1 A\ A AA re 

4iu; O.OOo 


1 TA 


|CCAA\ AA A77 

|oOO) 99. 9( 1 


1 A 1 A\ A AAfi 

|4iu; O.OOo 


o ceo 
2.0O 


lAAl \ Al A 

|001) 91. U 


1 Al A\ OA 

|oio; 8.0 


Z. (0 


l£JAA\ AA A'TO 

|dOO) 99.9/3 


1 K 1 A\ A AAO 

olo; O.OOo 


O OA 

3.29 


1 1 Al \ K O ^ 

|10i) 08. < 


1 Al A\ OO A 

uio; 22. u 


3.81 


101 fl\ 58 Q 


101) 21.7 


4.11 


|201) 43.8 


110) 18.8 


4.57 


|002) 64.0 


|011) 17.1 


4.71 


|110) 26.4 


|201) 16.9 


4.93 


|301) 30.6 


|201) 17.2 



Table IL Angular momentum basis content at a^c ~ 
0.36\/ft//itJ of the eigenstates of the full Hamiltonian 
for ^■^Rb and ^Li in a harmonic potential with f3 = 
(mLi<j-'ii)/(mR,bCJ^[,) = 0.4. The molecular states are almost 
pure product states with n = 0, L = and a well defined CM 
excitation A'' whereas for all other states the CM dynamics 
and the relative dynamics are entangled. 

small sinusoidal oscillations of the magnetic field around 
a fixed value B' close to a Feshbach resonance: 

Bit) ^ B' + bsm{LUBt). (34) 

If hujB matches the energy difference between two sta- 
tionary states, it is possible to transfer population from 
one state to another (for instance to the molecular state) 
or to create well controlled quantum mechanical super- 
position states. 

The modulation technique has also been used to pro- 
duce and make spectroscopy on ultracold heteronuclear 
®^Rb-*^Rb molecules to probe the excitation spec- 
trum of a Fermi gas in the BEC-BCS crossover regime 
[Slj . and to infer the lifetime of short-lived Feshbach 
molecules by looking at the width of the resonance 34] . It 
differs in princi ple from the traditional RF spectroscopy 
used in e.g. [tI. IBSI. [36j in that the atoms arc not trans- 
ferred to other hyperfine levels and in that it is the tun- 
able interatomic interaction that drives the transition. 

When the magnetic field is varied in time, the inter- 
action Hamiltonian is explicitly time dependent, and we 
have different options to determine the evolution of the 
system. 



A. Weak amplitude modulations of B-field 

If the B-field modulation is weak, it constitutes only 
a small time-dependent perturbation on the system, and 
a calculation naturally departs from an expansion in the 
adiabatic eigenstate basis. We label the complete set of 
eigenstates of the Hamiltonian at a given field strength 
by \'^q{B)) and their energies Eg, and we expand the 
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wave function 



(35) 



Inserting into the time-dependent Schrodinger equation, 
and projecting onto the orthogonal states {'ip{B{t))\ we 
get the coupled amplitude equations 



dt 



c,it). 



(36) 

If we take the sinusoidal magnetic field dependence ([34]) 
in we get the following equation for the coefficients 



ih 



dcp(t) 
dt 



Ep{t)Cp{t) - ihC0s{uJBt) ^ VLpq{B{t)) Cq{t). 

(37) 



where the Rabi frequencies U,pq are given by 



d 



^pq{B) = buB ^(*p(S)|*,(S')> 



(38) 



B'=B 



For a weakly modulated field, & <C A, the time de- 
pendence of Ep and Q.pq can be neglected and the mod- 
ulation is equivalent to, e.g., the radiative coupling of 
atomic energy levels, li fiWB is close to the energy differ- 
ence between two eigenstates p and g, and we have all the 
population in the states p initially, we therefore expect 
to see Rabi oscillations between the two states 



|c,(t)r = i-Mor=^g^ 

where 



sm^ f'^^l (39) 



n'p^ = ^nl^ + {ujB-{Ep-Eq)/hf. (40) 



B. General dynamical equations in modulated 
B-field 

The expansion ([55]) is general and may be applied for 
any modulation of the magnetic field, but since the eigen- 
states in that expansion have to be found first by numer- 
ical diagonalization of the coupling Hamiltonian, we find 
it more convenient to use instead the basis (|12p . where 
only the relative coordinate part of the wave function 
depends explicitly on time through the scattering length. 
We then get 



ih 



dcp{t) 
dt 



(41) 



At) 

(42) 



where Hpq is given in p9|) . 

The explicit expression for the time derivative of 
{ti}p{t)\ipq{t')) in the basis ^ is 

^^{^p{t)\^q(t')) 

= ^ {LpMplpmp\Jmj){LqMqlqmq\Jmj) 



/ 



(43) 



where we have exploited the fact that only the relative 
coordinate wave function ip is time dependent. We note 
that our use of the adiabatic eigenfunctions of iJrci re- 
moves the need to numerically compute any matrix ele- 
ments of the singular interaction potential, and thus we 
make optimum use of our knowledge of the analytical, 
irregular solutions to the pseudopotential equation. 

It is a great advantage of the pseudopotential model 
that the matrix element {(pi,oo\'fiu'oo) can be calculated 
analytically, and we can hence express the time derivative 
using the chain rule: 



d 



^(</?np00(i)IV'n,00(i')>r 



dB da^ 



-7(</5<.pOo|'Pi/'oo) 



where dagc/dB is given by 

da.r -A 



dB [B - Bo)2 
and dasc/dv is given by 

dflsc _ Osc _ Obg B ~ Bo — A 



with 



diy [/(i.)]2 [f{,y)Y 



B^Bo 



(44) 



(45) 



(46) 



(47) 



The overlap ('Pi/oo|¥'iy'oo) is provided in the Appendix and 
all coefficients in our coupled set of equations are thus 
given by analytical expressions. 



Results 



According to the discussion in Sec. IIV Al a small am- 
plitude harmonic oscillation of the magnetic field at a 
frequency such that huj is resonant with the energy dif- 
ference between two eigenstates will inevitably lead to 
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full contrast Rabi oscillations as long as there is a finite 
coupling between the two states. For J = such a finite 
coupling arises whenever the basis expansion ((33|) of the 
two eigenstates at the given magnetic field B' contains 
basis functions of the form <i>Aroo(R)<Piyoo(r) respectively 
$Aroo(R)¥'i.'oo(r), see Eq. HI 

This means e.g. that efficient association of heteronu- 
clear molecules should be possible. We can illustrate it in 
principle using the ^^Rb-^Li system treated above. The 
scattering length at zero magnetic field was measured by 
C. Silber et al. ^ to be 20(+9/-6) Gq. A Feshbach reso- 
nance with this background scattering length gives rise to 
the energy spectrum shown in Fig. [2l If we oscillate the 
magnetic field resonantly around B = Bq — 0.05A (corre- 
sponding to flsc/aiio = 0.36) we can produce full contrast 
Rabi oscillations between the lowest non-molecular state 
and the molecular ground state (illustrated with the low- 
ermost arrow in the spectrum). The dynamics is shown 
in Fig. [3] In practise one might want to use an adiabatic 
pa ssag e to make a robust full transfer as demonstrated 
in [Til . 



3 
UJ 
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0.15 



Figure 2: (Color online) Spectrum for *^Rb-^Li in an op- 
tical lattice well with u-Rb = 2-n x 200 kHz and jS = 
mLitt'Li/("^Rb'i'Rb) = 0.4 in the vicinity of a Feshbach reso- 
nance with flbg = 20 oq. Molecular states that have a CM 
excitation (A'^ > 0) are not shown. Association of a heteronu- 
clear molecule is indicated with the lower, red arrow while 
transfer to a state with angular momentum excitation is il- 
lustrated with the upper, green arrow. 



Another interesting fact is that the coupling of the 
CM and relative motion makes it possible to populate 
states with nonzero CM and relative angular momentum 
by oscillating the magnetic field close to a Feshbach res- 
onance. It is intriguing that in this way a purely central 
interatomic interaction can be used to put angular mo- 
mentum into the relative motion at the cost of exciting 
the CM motion as well. 

For example the uppermost arrow in Fig. [2]illustrates a 
transfer from the non-molecular atomic ground state to a 
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Figure 3: (Color online) Resonant association of heteronuclear 
molecules by an oscillating magnetic field. The time depen- 
dent magnetic field is B{t) ^ Bo - 0.05A -I- 0.005A sin(i.^st) 
where u>b = 2ii x 2660 kHz. The blue line shows the popula- 
tion of the atomic ground state whereas the population of the 
molecular ground state is shown with the red, dashed line. 
LjRb = 27r X 200 kHz, (3 = mLitj£i/(mRbijib) = 0.4. We have 
included states up to A'' = 3, n = 3 and L — I — 6. 



state with 59 % content of L = / = 1 (see Tab. |n|. If the 
magnetic field is oscillated resonantly, Rabi oscillations 
between these two states can be produced (Fig. HJ and 
again a 7r-pulse or an adiabatic passage could in principle 
be used to make a full transfer. 



V. DISCUSSION 

In summary, we have found the eigenstates of motion 
of a pair of harmonically trapped, different atoms inter- 
acting via a short range central potential. We have inves- 
tigated the dynamics of the system when the interaction 
strength is modulated, by changing the strength of an ap- 
plied magnetic field, and in particular we have found that 
resonant transitions between the molecular and atomic 
states can be coherently driven in this system. Interest- 
ingly we find that even though the interaction is central, 
states with non-vanishing CM and relative angular mo- 
mentum of the atoms can be produced. Although our 
calculations were done explicitly for the pseudo-potential 
P^ . this is not a formal restriction on the separation 
of the problem and the methods applied, and other cen- 
tral interaction potentials are readily treated by a similar 
analysis. 

We believe that the transition dynamics analyzed in 
the present paper may be studied with current experi- 
mental techniques. Our work identifies a way to produce 
molecules without crossing the Feshbach resonance, and 
we imagine that the states with nonvanishing relative 
angular momentum may lead to interesting momentum 
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Figure 4: (Color online) Rabi oscillations between the atomic 
ground state and a state with 59 % content of L = / = 1. 
The time dependent magnetic field is B{t) — Bo — 0.05A + 
0.005Asm{ujBt) where ujb = 2n x 601 kHz. The blue line 
shows the population of the atomic ground state whereas the 
population of the excited angular momentum state is shown 
with the red, dashed line. The other parameters are as in Fig. 

m 



A. Radial integral for harmonic oscillator wave 
functions 



The normalization constant of the harmonic oscillator 
wave functions is given by 



1 



1V2 



_ r(n + ; + 3/2) 
~ 2n! 

where we have used ([H 22.2.12, jilj 7.414.3): 



C r(a+«+i) ^ = „Rea>0 
I 0, m 7^ n, Re a > — 1 



(50) 



(51) 



distributions [H, |3^ |4^ of the individual species when 
the atoms are released from the lattice. 
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VI. APPENDIX: CALCULATION OF MATRIX 
ELEMENTS 

This appendix derives the matrix elements used in the 
Feshbach molecule problem. We consider the radial wave 
functions for an isotropic 3-dimensional harmonic oscil- 
lator (dlD 

i?„,(r) = 7V„iLi+i/2(r2) r'e-'^'/2 ^^g) 

where L is a generalized Laguerre polynomial, and the s- 
wave solutions to the problem of a particle in a harmonic 
oscillator modified by a regularized s-wave J-function po- 
tential at the origin ([?f|) 

i?.(r) = — 2-'^-i/2i^2.+i(2), z = V2r (49) 
r 

where D is the parabolic cylinder function. N^i and N^, 
are normalization constants. 



The radial part of the dipole matrix element between 
oscillator states involves the integral 

/■OO 

= i / 4+i/2(u) Lj,+i/2(u) u(''+'+2)/2e-"du 
2 Jo 

(52) 

Since the angular integral has the selection rule \ td\ = 1 
we can, without loss of generality, assume that I' = I + 1 
and employ the recursion relation 

LZ-\u) = LZiu) - LZ_,{u) (53) 

for the generalized Laguerre polynomials (^^ 22.7.30) to 
obtain 

L:/^/^(«) L5+^/^(.) .(''+'+2)/2e-"du 

irv+l+5/2) 

"2 (^\ K^r,' ,n - 5n' ,n-l) (54) 

where (jSip was used in the last step. Including the nor- 
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malization ([50)1 we obtain 

Rn' ,i+i{r)rRni{r)r^di: 



l T{n' + 1 + 5/2) nl 
r(n + ; + 3/2) (n')! ^ "''""'^ 



= \/n + l + 3/2 Sn'^n - Vn 5n',n-l (55) 

where the identity r(a; + 1) = a; r(a;) has been used. 

B. Radial overlap of parabolic cylinder functions 

The normahzation constant of the parabohc cyhnder 
wave functions (P^)) is given by 

1 /"OO /"OO J 

^(_^)_^(_^_l/2) 



■,-2i/-3 



r(-2i^-i) 

TT v(-i^) - V'(-'^ - 1/2) 
2 r(-^- i/2)r(-i/) 



(56) 



where [4l| 7.711 has been used to solve the integral and 
the duplication formula for gamma functions |29j 6.1.18 
has been used in the last equality. ■0 is the digamma 
function: 



(57) 



In the case i/' ^ v we apply [4l| 7.711 to determine the 
overlap integral 



2{v' - v) 



1 



1 



T{-v')T{-v - 1/2) T{-v)T{-v' - 1/2) 



^r(-i.' - i)r(-^') ^T{-u-\)T{-i,) 



(58) 



where 



(59) 



C. Dipole matrix element between parabolic 
cylinder and harmonic oscillator p-wave functions 

We evaluate 



R,{r)Rni{rydY 



r'^dr, 



(60) 
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The generalized Laguerre polynomials are given by ([2^ with v = —{2v + 1) to get 
22.3.9) 

K{u)=Y.-'^u\ a, = {-lf[l^_l^l^^ (61) i?.(r)i?„iMrMr 

and we integrate term by term using ([H 7.722 for u > - /v /V , ^ V 2^2fc/ i \k ( + ol \ }_ ^{'^k + A) 
0): ^16^0 ^ H"-^ /fc! r(fc + 2-;.)- 

(63) 



^ ' r(/i/2 + ;>/2 + 1/2) ^ ^ 
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